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The generation of harmonics of the voltage response is considered when an AC current is applied 
through a superconducting film above T c . It is shown that almost at all temperatures the mechanism 
of the temperature oscillations created by the AC current and the temperature dependence of the 
resistance dominates over the isothermal nonlinear electric conductivity. Only in a narrow critical 
region close to T c the latter is essential for the generation of the harmonics. A detailed investigation 
of harmonics generation provides an accurate method for measuring the thermal boundary conduc- 
tance between the film and the insulating substrate. The critical behaviour of the third harmonic 
will give a new method for the determination of the lifetime of metastable Cooper pairs above T c . 
The comparison of the calculated fifth harmonics of the voltage with the experiment is proposed as 
an important test for the applicability of the employed theoretical models. 



I. INTRODUCTION 

Third harmonic generation is frequently used in a lot of scientific fields to obtain information which is impossible to 
extract fcpm the total signal or even from the first harmonic. Third harmonic techniques are particularly developed 
in opticsfil, and in the physics of semiconductorim. It is also possible to explore superconductor properties with such 
techniques as was already done for thermal properties in bulk and thin film materialsuu, for magnetic responacEl, 
non-linear effects in YBa2Cu307_,5 bicrystal Josephson junctional and the nonlinear response in the vortex stateo. 

The aim of the present paper is to give a theory for the harmonics generation of the voltage response when an 
AC current is applied through a conducting film deposited on an insulating substrate. Systematizing the results for 
the temperature oscillation mechanisnoQ we describe how the effects of isothermal nonlinear current response can 
be extracted. Experimental data processing of the harmonic response can give not only an accurate method for the 
determination of the thermal boundary conductance G(T) but also the lifetime constant To of metastable Cooper 
pairs above T c . We suppose that the geometry consists of a strip of a superconducting film. 

Generation of 3rd and 5th harmonics in the voltage response of such a sample due to thermal effects of the electric 
current is presented in Sec. Q. In Sec. HI is performed the comparison of the thermal mechanism with the isothermal 



mechanism of nonlinear electric conductivity^ and it is described how the effects of these two mechanisms can be 
separated. Finally in Sec. ^ it is considered how the investigation of the 3rd and 5th harmonics can be used for 
practical applications in the physics of bolometers and for the determination of the parameters of superconductors 
important for understanding the fundamental processes in these materials. 

II. MODEL 

We consider the temperature T of a strip patterned from a thin superconducting film with length L and width w, 
deposited on an insulating substrate when a bias current I is applied. Due to ohmic heating, the film temperature T 
is slightly higher than the substrate temperature T su b: 
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T(t) = T sub + AT(t). 



(2.1) 



We suppose that the current is low enough to ensure very small temperature variations A T -C T. In our model we 
will use the phenomenology of the thermal boundary resistance G(T) between film and substrate in which the total 
power of emitted phonons is defined as: 



P = GAT: 



(2.2) 



for a general review of the thermal boundary resistance between solids see for example the Ref.0 The power P 
is proportional to the temperature increase AT and to the total thermal conductance between film and substrate 
G = wLg(T), which is equal to the thermal conductivity g(T) multiplied by the cross-section wL of heat exchange. 
We suppose that the temperature in the film is almost constant across the film in all directions. Taking into account 
the thickness of the film dgi m and thermal conductivity k of the superconductor, we assume that the thermal boundary 
resistance is much higher than the thermal resistance of the thin film (dfii m /K <C l/.<?)- Due to the insulating nature 
of the substrate the interface heat exchange is a purely phononic phenomenon and we expect the thermal boundary 
conductivity g(T) to be a smooth function of the temperature, and also G(T) w G(T su b). We assume that the 
frequency of the AC current is low enough, u> -C G(T)/C(T), where C(T) is the heat capacity of the film. In this case 
effects related to the specific heat are negligible and the temperature follows instantaneously the current through the 
film, so that T = T(I). In these conditions, the total power is expressed by Ohm's law 



P = R(T)I 2 , R{T) 



P(T)L 



/ = wd 



film J 7 



(2.3) 



where p and j are respectively the resistivity and the current density. The comparison with Eq. (| 
power gives the temperature difference 



for the emitted 



AT = 



R( T ) _j2 = P( T ) d fi\m , 2 
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g(T S ub) 



(2.4) 



Taking into account that the resistance of the film is temperature dependent R(T) and using Eq. (2^) and Eq. (p.l|), 
we easily derive the power expansion of the resistance as a function of the current: 



R(T) = R(T suh + AT) = i? T sub 
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where 



R'(T) = ^iZ(T), R"(T) = ^R{T). 



(2.6) 



Eqs. (2.2-2.5) are essentialJbr the physics of thermo-electric oscillations in HTSC fi m structures0 and the d agnostics 
of HTSC microbolometersEj. Then for the voltage, we have: 



U = R(T)I w RI 
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Let us now consider the case of a harmonic current 



I{t) = I COSUjt, j Q = Io/wdfiln 



Using trigonometric relations 



3 1 o 3 

COS bjt — — COS ibjt + — COS LUt 

4 4 

cos 5 ojt = — cos Suit + — cos 3u>t + — cos cot 
16 16 16 



(2.9) 



and writing the voltage in a Fourier series 
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U (t) = U\ u cos uot - 
we obtain for the amplitude of the harmonics: 
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Analogous formulae fen-the Zlo method applied to thermal conductivity measurements were reported for rod- and 
filament-like specimens Ji-3 For the case of linear dependence of the resistance vs temperature well aJjflve T c for some 
layered cuprates R(T) = Ui u /Iq « AT + B, these expressions up to Iq -terms were analyzed in Refs.CfB in the context 
of the determination of the thermal boundary conductance G{T) of the YBaaCuaC^-^/SrTiOa interface explored by 
the third harmo nic ge neration of the voltage across the strip. In this case cf.Qa R' = A, R" = and for small enough 
currents Iq <C \J G/A we have the approximations: 
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For experimental data processing of temperature sweep measurements we can express all measurements as functions 
of the time t: 



R\T) = | 
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(2.13) 



where a dot means time differentiation; this is helpful if we use numerical filtration of the signal. We suppose 
that for small enough currents Iq —> 0, the intensity of the harmonics is decreasing fast with the overtone index 
\Uiu\ S> \Usu;\ 3> |Cso,[. In this case at fixed temperature we can perform a polynomial fit of U nu / 'RIq versus Iq, and 
use at most a second-degree polynomial for n = 1, 3, 5. 

If we investigate the temperature dependence, roughly speaking U\ u determines the resistance R(T), and, subse- 
quently, Uj, u the thermal conductance G(T). The comparison of the measured with the prediction of Eq. (2.11) 
can be considered as an important test for the consistency of the model. Using simultaneously the recordings of the 
first, third and fifth harmonics of the voltage we can obtain better accuracy for the experimental determination of 
the resistance and the thermal conductance as follows, 
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where E nu = U nuJ /L are the amplitudes of the electric fields for n — 1, 3, 5. These high precision methods could be 
helpful for extracting the electric field nonlinearity analyzed in the next section. 



III. NONLINEAR CONDUCTIVITY: NUMERICAL EXAMPLE 



Even for negligible temperature oscillations the nonlinear current response of the fluctuation conductivity can 
generate harmonics. For superconductors the coefficient A in the nonlinear term of the current 

J = ^--A(T)E3 + 0(E 5 ) 7 E = p[j + (pjfA(T)]+0(f) (3.1) 

describes the depairing effect of the electric field. In other words, an electric field accelerates the metastable Cooper 
pairs above T c and this acceleration increases the decay rate and finally the density of fluctuation Cooper pairsli3. 
The nonlinear conductivity coefficient A has a strong critical singularity as a function of the reduced temperature 

T — T 

e=i_ii:«l. (3.2) 
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For a thin film, for example^, we have 



A(e) 
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(3.3) 



where e is the electron charge, T c the critical temperature, £(0) the coherence length and r re i is the ratio of the lifetime 
constant tq of the fluctuation Cooper pairs to the theoretical value 
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(3.4) 



derived within microscopic theo ry in low-coupling and negligible-depairing approximation. 

According to Ohm's law Eq. (2.3), and Eqs. (2.9) and (2.11) both considered mechanisms of 3rd harmonics gener- 
ation take the form: 
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(3.5) 



where p'(T) = dp(T)/dT and the superscript denotes the thermal (T) or electric-field (E) origin of the harmonics 
generation. is the total voltage signal measured by the lock-in. It consists of two contributions: the thermal part 
is proportional to the thickness of the film because heat is generated in the bulk of the film but leaves it through 
the interface; the electric part is thickness-independent. The former behaves near T c as a constant with a correction 
proportional to e -2 , while the latter displays a strong singularity proportional to e~ 4 , as given in Eq. (3.3). 

Let us make an order of magnitude evaluation taking the variable T to T c and taking to as given by Eq. (|3.4|), i.e. 
Trei = 1: 
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where pn(T c ) is the extrapolated normal conductivity from temperatures well above T c . For high-T c cupratesSQ 
have 



we 



g(T c ) « 1000 W/Kcm 2 , £(0) = 1.1 nm, pjv(lOOK) = 132 /xficm, 
p^r(lOOK) = 1.09 /xficm/K, T c = 90 K, d film = 50 nm, e ET ^ 1%. 



(3.7) 



This means that the nonlinear effect due to the electric field will dominate in a 1 K region around T c , and will 
still be observable 2 or 3 K above T c . In order to extract this term we have to determine G(T) from Eq. ( |2.14 ), 
using the experimental data in the temperature interval, e.g. , (1.1 T c , 2T C ), well above the critical region. Being a 
purely phononic process for insulating substrates the thermal conductivity g(T) has no critical singularities and can 
be reliably extrapolated to T c using the polynomial fit of the experimental data. Then this polynomial approximation 
Gfit(r) = Lwght(T) can be used in Eq. (2.11) together with the data for the resistance in order to compare the 
experimental data with the prediction from the model of the thermal oscillations. The model has two functions 
g(T) and p(T) which have to be experimentally determined and fitted. We suppose that the simple model of purely 
thermal oscillations will describe the main qualitative property of the temperature-dependence of the 3rd harmonic 
- the maximum at T c determined by p' . This maximum already contains a critical behaviour due to the critical 
behaviour of the conductivity l/p(T) which has a significant fluctuation part close to T c . The singularity created by 
the depairing effects, however, is stronger close to the critical region and we can determine this electric nonlinearity 
using the experimental data and the fitted thermal conductance. Taking into account Eqs. (2.11) and ( |3.l[ ) we derive 
the experimental definition of the nonlinear voltage and its coefficient 
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where 1^'^ is the small parasite amplitude of the 3rd harmonic of the current created by the current generator. The 
current I should be small enough 
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(3.9) 

in order for the Jg term to be much smaller than the ig term, and for the 7q/GL correction to be negligible. The 
nonlinear coefficient determined in this way should be compared with the results of the fluctuation theory for a layered 
superconductor^ 

3 i 3„,2 i , 5 ^31 
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s6(0) [e(e + r)] 7/2 

where £ o (0), £b(0) and £ c (0) are the coherence lengths, s is the interlayer distance, and r = (2£ c (0)/s) 2 . We suppose 
that the electric field is applied along the a-axis. For high-T c cuprates a significant nonlinear effect can be observed 
only for very thin high-quality films. . .. . 

Very strong nonlinear I — V characteristics V oc I 3 are observedOlla at Kosterlitz-Thouless transition temperature 
Tkt- At this temperature 3rd harmonic will have a A-shaped maximum analogous to the maximum observed at 1st 
harmonica. For temperatures slightly above Tkt the 3rd harmonic is created mainly by fluctuation Cooper pairs. 
Except in the considered narrow critical region the thermal oscillation mechanism will give the main part of the 
harmonics generation in the general case. Only for some conventional superconductor films the cet parameter may 
be large enough so that the thermal contribution can be neglected but this problem requires additional investigation. 

IV. DISCUSSION AND CONCLUSIONS 

Looking at the expression of €et, it appears that in high-T c superconductors the electric field nonlinearity for the 
generation of the third harmonic can be with better accuracy measured for high-quality films as thin as possible. 
In this case the investigation of the third harmonic between T c and T a few degrees above T c can be used for the 
determination of the lifetime constant To of fluctuation Cooper pairs. 

Far above the critical region, we considered a quite universal mechanism of the voltage harmonics generation due to 
thermal oscillations. This mechanism is very effective close to the phase transition, superconducting or metal-insulator, 
where the resistivity changes significantly in a narrow temperature interval. Perovskite thin films, high-T c cuprates 
and manganates with colossal magnetic resistance are very promising for technical applications like bolometryD'B, 
for example. Detailed investigation of thermal properties by simple measurement of 3rd harmonic can give useful 
information about the quality of the sample. A detailed investigation of the temperature dependence of harmonic 
generation can give a precise method for the determmatiorOJ of the thermal boundary conductivity <?(T). Very 
intense third harmonic generation could be expected for F^S^CaC^Os whiskers placed in vacuum. The black body 
radiation is analogous to the phonon radiation in the insulator but the emitted power is much smaller, because at low 
temperatures the power is inversely proportional to the square of the wave velocity in the heat sink. 

In such a way the fundamental problem of the physics of superconductivity, the determination of the Cooper 
pair lifetime constant and the problem of the determination of the thermal boundary conductivity <?(T), important 
for many technical applications, can be solved simultaneously by investigating the voltage harmonics generation in 
superconducting microbridges. 
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APPENDIX 



In this appendix we will consider in short some formulae, convenient for programming, for smoothing and numerical 
differentiation of a function measured experimentally fj = f(tj) for j = 1, 2, 3, . . . , N with some error. The main 
idea is to perform in the beginning the interpolation of the function using the Taylor expansion at some point t and 
the values of the function averaged with respect to the experimental data and its 1st and 2nd derivatives (/(i)}, (f{t)) 
and </(*)} 



fit) « (/(*)> + </(*)>(*-*) + 2</(*)>(* " *) 



(4.1) 



Applying this parabolic interpolation to the interpolation points t — U we will determine the averaged values of the 
function searching the minimum of the sum of squares of the errors 
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where the weights Wk > depend on the distance from the interpolation point tk to the center of the interpolation 
t, and for simplicity we skip this argument in the notation for the averaged function (/) = (f(t)), . . .. One possible 
weight function is 
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where W T is the half- width of the interpolation interval and S T is the sharpness of this window. The minimum of the 
Gaussian function G with respect of (/}, (/), (/) gives the equations 
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which have the solution 
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The formulae simplify significantly if we use equidistant points tj = jAt and take the center of the interpolation at 
one of them t = U, i.e. 



.n+k^(h) + (.n)Atk + -{.mAt) 2 k 2 . 

In this case the odd momenta are zero, [t] = and [r 3 ] = 0, and we have 
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A further simplification is to take a sharp window for interpolation 
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This means that sums include only 2L + 1 points centered at i and we can use directly the integer index k instead of 
the argument r 
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The simplest possible example is probably the 7-point averaging for L = 3, in this case for 2 < i < N — 2 Eq. ( |4.5| ) 
gives: 
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[3(/i+a - fi-s) + 2(/ i+2 - /i-a) + (/i+i - /i-i)] , 
[5(/i- 3 + /i+a) - 3(/i_ a + /n-a) + 4/,] . 



(4.13) 



Some of these formulae can be used sequentially for smoothi ng, d ifferentiation and further smoothing of the derivative. 
For the ends we can apply the parabolic interpolation Eq. (4.8) 



</*} = (ft) + (fi)At(k -i) + -(.h)(At) 2 (k - if, 



(4.14) 



for i = 3 and k — 1, 2. The same formula has to be appl ied for the right end of the data for i — N — 3 and 
k = N — 1, N — 2. Finally the parabolic interpolation Eq. (4.1) can be used for tabulation of the inverse function t(f) 



(fit)) 
</(*)> 



l 2(f(t))(f-(f(t))) 1 
if(t)) 2 



(4.15) 
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